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Abstract 

We discuss self-similar solutions to 0(4) textures in Minkowski 
space and in flat Friedmann-Robertson- Walker backgrounds. We 
show that in the Minkowski case there exist no solutions with wind- 
ing number greater than unity. However, we find besides the known 
solution with unit winding number also previously unknown solutions 
corresponding to winding number less than one. The applicability of 
the non-linear sigma model approximation is discussed. We point out 
that no spherically symmetric exactly self-similar solutions exist for 
radiation or matter dominated FRW cosmologies, but we find a way 
to relax the assumptions of self-similarity that give us approximative 
solutions valid on intermediate scales. 
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1 Introduction 



Since some time, it has been realized |l|] that defects (textures) associated 
with the non-trivial winding of massless scalar fields may be of interest even 
though they are intrinsically unstable if the winding number is large enough 
1^, 1^, ^. Indeed, it is the fact that textures continually enter the horizon 
during the evolution of the universe that makes the spectrum of density 
fluctuations near scale-invariant (although not Gaussian) and makes textures 
promising candidates for large scale structure formation even in light of the 
COBE resuhs |, i- 

A simple theory that admits global textures is given by the lagrangian 

where $ is a 4 component real scalar field. Here 0o is the sjTumetry breaking 
scale and A is a dimensionless coupling constant. For convenience we do the 
rescaling ^ ^ (pQ^, so the action for this theory will be 

Sm =<PlJ d'x^g {\d,^ ■ - f ($^ - If) , (2) 

were w = A0q and g is the determinant of the space-time metric. 

Upon quantization we have in this theory one massive Higgs particle with 
mass tuh = \Jw/2 and — 1 massless (Goldstone) bosons. We will however 
not be interested in the particle spectrum, instead we will only consider the 
classical equations of motion, for the field 

d,{V^d^^) = -^gw{^^ - 1)$. (3) 



For cosmological applications, it is important that the Goldstone modes 
remain massless, creating long-range correlations and field dynamics over 
cosmologically relevant length scales. Arguments, based on quantum gravity 
effects, have been given |^ which seem to make the survival of exact global 
symmetries questionable. This statement is, however, based on unknown 
physics at the Planck scale. More specifically, mechanisms present, e.g. in 
string models may well protect the potential of the Goldstone modes of the 



texture scalar fields. (For a recent discussion of such mechanisms, see |T^ .) 

Here we will not enter into this discussion but simply assume that tex- 
tures can exist and study the properties of their evolution when the dynamics 
is given by the action (0). 

We will parameterize the field using hyper-spherical coordinates p, X; ^ 
and (p, in the following way, 

(p{r,t) = p(cos X, sin X cos 6', sin x sin 9 cos sin x sin 9 sin (p) . (4) 
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We will look at the "spherically symmetric" (or hedgehog) ansatz, were 
we let the coordinate functions depend on time t and the spatial spherical 
coordinates r, 9 and as p = p(r, t), x = x(^) t), 9 = 9 and (p = (p. 

It is common to consider the field as a stiff source, which means that one 
is assuming that the self-coupling of the field is much stronger than the self- 
gravitational coupling. Thus only the background metric is required in the 
equation of motion for the field. The perturbation of the background metric 
can then be calculated from Einstein's equations with the stress-energy ten- 
sor of the texture field added. The applicability of the stiff approximation 
in the self-similar case is discussed in . 

When one studies the formation of large scale structure in the early 
universe, the background metric is taken to be Friedmann-Robertson- Walker 
(FRW). We will discuss how one can find solutions valid at medium large 
scales using the equation for the Minkowski background. 



2 Minkowski background 

For a Minkowski background the equations of motion (§) in terms of the 
hyper-spherical coordinates become 

r^ip'^x)' = (p'^r'^x'y - sin 2x , . 

(5) 

r^w{p^ — I) = r'^ix^ — x'"^) — 2 sin^ x + ((^^P')' ^ ^^P)/P- 

For length scales larger than the inverse mass of the radial "Higgs" mode 
= (A0g)^^/^ the dynamics of the field can be described by a nonlinear a 
model (NLSM) [0. The NLSM is characterized by that the field is exactly on 
the vacuum-manifold everywhere, thus p = 1. In this case the first equation 
of admits a self-similar ansatz y = r/t and becomes 

(1 - y^){y\yy + 2yxy) = sin2x(l/). (6) 

This equation has a singular behavior at ?/ = and y = ±1, the conditions 
for regular solutions are sin2x(0) = and sin2x(±l) = 0. 

This equation has the well known solutions found by Turok and Spergel 

i- 

x{y) = rrni ± 2 arctan(±?/), (7) 

where m is an integer. 

These solutions are indeed very special, coming from the spherically sym- 
metric self-similar ansatz to the non-linear sigma model approximation. An 
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important feature is, however, that these symmetric scahng solutions ap- 
pears in the numerical simulations as attractors 0. In fact, many of the 
calculations of the effects of textures, e.g. on the microwave background 
radiation, rely on the use of these simple analytic solutions |1^. 

In order to see when the NLSM approximation is applicable we go a 
step beyond it. Instead of insisting on p = 1, we will assume only that 
the derivatives of p are negligible in the equations of motion. We will then 
recover the NLSM equation for %, so for x we will use the selfsimilar NLSM 
solution (|^. The second of the eqs. (||) p becomes 



p" = 1 + {X^ -X'^ -^sm^x)/w, (8) 



which gives 

upon insertion of our self-similar NLSM solution for x- K can be checked 
that p' and p can be neglected in the first of the eqs. (|^) if + >> 1/w. 
Thus we conclude that there exists a ro >> such that the selfsimilar 

solution (|^ is vahd for all + > tq. 

If we have a field that initially for t<0, i(:^>>l/wis described by 
x(r, t) = 2 arctan(— r/t) we will have an unwinding event for + < Tq 
where the field is forced to leave the vacuum manifold. The solution (0) is 
valid right to the time t = r — ro, when the information from the unwinding 
event reaches r. We thus can match the solution at t = 0, xl'"? 0) = vr, 
x(r, 0) = 2/r, with = 27r — 2 arctan(r/t), valid for < t < r — Tq. 

How the field behaves for t > t — tq depends on the details of the unwind- 
ing event and must be decided by making a numerical simulation of the full 
field equations One thus find that the field after the unwinding goes 



asymptotically to the NLSM solution xl'"?^) = vr + 2 arctan(r/t) for t > r. 
This solution describes an expanding shell of goldstone bosons, the winding 
number is zero. 

Before we continue discussing Minkowskian self-similar solutions, we study 
the equations for the FRW background metric. 

3 FRW background 

We will discuss the evolution of spherical textures in a fiat FRW background 
ds^ = a\r]){dr]^ - dr^ - r^dO^ + sin'Od^^)). (10) 
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Here rj is the conformal time. The time dependence for a is a{ri) oc 77", where 
a = 1 corresponds to a radiation dominated universe and a = 2 corresponds 
to a matter dominated universe. 

The equation of motion for the NLSM is now 

('^^Xr)r - ?^^(x,,r, + 2^Xr,) = sin 2x (H) 

We are interested in solutions where unwinding can occur so we try the 
selfsimilar ansatz xi^^t) = xl'"/^)? where t = r/ — 77^, and rj^: is the time 
for the unwinding. The equations with a 7^ admit this ansatz only if 
?7* = which coincides with the time for the big bang singularity. One 
would nevertheless hope to have some use of this ansatz if we are interested 
in expanding textures that unwinded very early. 
The self-similar ansatz y = r/rj gives 

y\l - y'')xyy + 2y(l + y\a - l))xy = sin 2xiy). (12) 

However, we will show that with a = 1 or 2 there does not exist any non- 



trivial solutions to (|l^) passing through y = 1 with finite derivative. 

To show this we use the regularity conditions that we get by letting y ^ I 
in (|T2]) and the first and second derivative of that equation. We thus have 
at y = 1 

2axy = sin 2x, 

(a - l)xyy + (3a - 2 - cos 2x)Xy = 0, (13) 
(a - 2)xyyy + (6^ - 9 - cos2x)Xyy + (6(a - 1) + 2cos2x)Xy = 0. 



For a = 1 we find (cos 2^(1) — l)xy{l) = 0. This gives Xy{^) = since 
cos2x(l) = 1 implies Xj/(1) = 0. For a = 2 we find 

(cos2 2x(l) - y cos2x(l) + y)Xy(l) = 0, 

the two solutions for cos 2^(1) are cos 2^(1) = so the only real solution 
is also here Xy(l) = 0. 

We see that in both cases for regular solutions we must have Xy{^) = 0. 
Since xiv) = n7i/2, integer n are solutions to (|T2D with Xyi^) = and 
sin2x(l) = we conclude that there does not exist any non-trivial regular 
solutions. 
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Instead we go over to discuss the validity of using the equation for the 
Minkowski background as the hmiting case when we look at small scales. If 
we substitute rj = rj^.. + 1 in ([TT|) and assume that t << r/* we get 

{r^Xr)r - r^iXtt + 2—Xt) = sin 2x. (14) 
V* 

Now if we can neglect the term linear in xt compared to Xtt we recover the 
Minkowski equation. Inserting the solution xi'^/t) = 2arctan(r/t) we see 
that this approximation is consistent only when |t|77* >> a(r^ + 1^). Thus 
we have to try a different approximation in order to get something valid for 
t ^ 0. 

We have found a way to get rid of the term linear in the time-derivative 
by a certain transformation. In the case a = 1 when the equation of motion 
is 

2 

ir^Xr)r - r^iXvv + -Xv) = sin 2x (15) 
we can make the substitution x(r, r^) = ^ilj{r,r]) and get 

(r'^ipr)r — '^^V'rjn = — sin2— (16) 
r]* V 

which become the Minkowskian equation after substituting rj = t]^ + t and 
neglecting t compared with 77*. Thus we find for a = 1 the solutions 

X{r,v) = -^■ipM{r/t) (17) 

valid for t << r]^^ where ipM{f/t) is any solution to the Minkowskian equa- 
tion. 

The same trick can be done in the case a = 2, but first we have to change 
to the coordinates u = Srjlr and t = rj^ in eq. (|ll]) giving 

{U\u)u - iT/nY/\^iXrr + 'Xr) = Sin 2^, (18) 

r 

which has the solutions 

x{u,r) = ^^JjM{u/t) (19) 
n + t 

valid for t << t^. 

In the linear approximation of Einstein's eqs. one can calculate the New- 
tonian gravitational acceleration from the unwinding texture solution (|^), 
the result is 
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where e = SttG^q, G the gravitational constant and f is the radial unit 
vector. This acceleration give rise to a velocity kick inwards of surrounding 
homogenous dust of the amount ne 0. 

For the solution (|^) the same calculation gives to first order in t/r*; 

9 = -e^,il-t/r.)r. (21) 

We notice that the acceleration ( pT]) is enhanced at t < compared with 
the ordinary (]20|), and vice versa for t > 0, which can be of importance for 
the form of the resulting matter perturbations. The velocity kick of the dust 
over the time interval {— to^^o}, to << is 25 arctanto/''", so at r << to it 
is still ne. 



4 New solutions 

We now examine the self-similar solutions of the Minkowskian equation of 
motion in greater detail. 

The solutions have a winding charge Q = ±1. One would perhaps 
believe that there exist solutions with higher \Q\ than 1 as has been claimed 
in P], but this is not the case. 

For a self-similar spherical texture we have |Q| < 1. This is implied 
by the theorem we prove in the appendix that if xiu) is a regular solution 
to (|^) with x(0) = then < \x{y)\ < ^r, for all finite y, and we have 
< |x(oo)| < vr. 

The argumentation of reference concerning solutions with |Q| > 1 is 
based on the erroneous assumption that there exist solutions satisfying the 
boundary conditions x(0) = 0, = n7r/2, with n > 1, (see the corollary 
of Lemma 1). 

We also want to emphasize that a boundary value problem such as (H) 
with x(0) = and x(l) = 7r/2, does not necessarily possess a unique solution. 

Actually we have by numerical means been able to demonstrate the exis- 
tence of whatseems to be a countably infinite set of additional solutions with 
total winding charge Q less than unity. These solutions are characterized 
by the number of oscillations around the value x = ^/2, and have rapidly 
increasing derivatives at the origin. 

We want to demonstrate the existence of these new solutions by accu- 
rate numerical techniques [|TT] with some modifications necessary to handle 



the singular points of the equation. We want to solve the boundary value 
problem using the shooting technique. The most straightforward strategy 
would then be to consider the initial value problem x(0) = 0, Xy{^) = Pi 
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and numerically integrate this to the point y = 1, we denote the solutions 
by xiUyP)- The equation xi^^P) = 7r/2 may now be solved for j3 by trial. 
However, because of the singularities it becomes numerically impossible to 
start the integration from ?/ = 0, so we have to modify our method. We 
must start the integration at a small distance away from y = and use a 
series expansion in order to get an appropriate initial condition. 

By making a series expansion of xiv) close to the origin of the form: 

oo 

X{y) = E a2k+iy"'^' (22) 

A;=0 

and inserting this into the differential equation Eq. (^), one can find the 
coefficients as, a^,... in terms of ai = p. One finds, e.g., 

2/3 - 4/3V3 
^ 10 

and 

3/3 - 3/33 + /35 

= 35 • 

Let us now consider the initial value problem, 

Xie) =Pe + asE^ + a^e\ Xy{e) = /3, (23) 

and integrate only up to y = 1 — e. For each e we may find a /3 such that 

x{l-e,P)+exy{l-e,P) = 7c/2, (24) 

we then must check that the value of /3 converges when we choose e smaller 
and smaller. It is also possible to use a more accurate extrapolation formula 
near y=l, one may again use the form of the original differential equation 
to write 

.1 ^ ^ 7 2 7 3, 7(1-7^) 4 , 

Y[l — z) = 7Z z z -\ z + ... 

^2^2 6 18 

where 7 = Xy{^)- This expression can also be used to continue the solution 
past the singular point y=l. 

With this careful treatment of the singular points of the differential equa- 
tion, its solution is otherwise straightforward. For the numerical solution we 
used a Runge-Kutta method with adaptive size control. We have employed 
this technique and thus discovered a set of such different /3's. In Fig. 1 the 
first four of the solutions corresponding to these /3's are displayed, we num- 
ber them with the mode number n starting with n = for the analytical 
solution. 
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These solutions appear to be very robust according to various stability 
checks we have made of our numerical algorithm, so we are confident in the 
belief that the presence of the solutions is not a numerical artefact. We have 
also checked that when we vary a in (0) around we still find solutions 
which approach our solutions in a continuous way when a ^ 0. From 
the conspicuously regular pattern of the first solutions shown in Fig. 1, we 
conjecture that the number of solutions is countably infinite. 

That self-similar solutions with winding number less than unity exist 
is potentially of great importance, since as shown in numerical simulations 
and backed by analytical arguments 0, configurations with Q > 1/2 
collapse and contribute to structure formation. We expect that well inside 
the horizon where spacetime is approximately Minkowski our new scaling 
solutions could play a dynamical role in structure formation. These solutions 
need a very high resolution numerical code to appear in the simulations 
since the derivatives at y = are very high. We plan to investigate these 
questions as well as the attractor nature of the solutions in future work. The 
applicability of the NLSM for the solutions with winding number less than 
one can be examined in the same way as for the analytical solution. The 
NLSM is valid for + > t\ and we find that we get a factor of around 
hundred extra in tq for each mode, the condition reads tq >> 100"/-^/?/;, 
where n is the mode number. We can use the solutions for times t <r — 
as for the analytical solution, but they can not be matched at t = r with 
any selfsimilar solution, thus for t > r the selfsimilarity will necessarily be 
lost. 

To conclude, we have investigated in quite some detail the nature and 
validity of the self-similar ansatz to the texture equations of motion. We 
have analyzed possible modifications when one goes beyond the non-linear 
sigma model approximation, the Minkowski background approximation, and 
the "ground state" arctangent solution. In future work, the effects caused 
by including the self-gravitational coupling will be investigated. 

We are grateful to P. Ernstrom for useful discussions. The work of L.B. 
was supported by the Swedish Natural Science Research Council (NFR) and 
EEC-SCIENCE contract no. SCl*-CT91-0650. 
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Appendix 



In this appendix, we prove the following theorem: 

Theorem If xiv) is a regular solution to with x(0) = then < 
\x{y)\ < 71") for all finite y, and we have < |x(c)o)| < vr. 
For the proof we need some lemmas: 

Lemma 1 If xiu) is a regular solution to (|^) with x(0) = then < 
\xiy)\ < TT for < y < 1. 

The proof is similar to one used in P] concerning the boundary conditions 
for the static ansatz xi^^t) = fir)'- 

We make the variable substitution x = l/y in (H) and get, 

(x' - l)xxx = sin 2x(x), x(oo) = 0. (25) 
We multiply this with Xx and integrate from x to cxd, 

- ^)xl 

Since Xx{x) = —^iXyi^) + C'(l/x)) for large x, ( PB| ) reduces to 

("OO 

- (a;^ - l)Xx(a;) - 2 / dxxxl = cos 2x(x) - 1. (27) 

Jx 

The left-hand side is always negative for 1 < a; < oo so we must have 

< \x{x)\ < TT for 1 < X < CXD. 

Corollary If x{y) is a regular solution to (^) with x(0) = and Xyi^) > 
then x(l) = vr/2, if x,(0) < then x(l) = -vr/2. 

This follows immediately from lemma 1 and the regularity condition 
x(l) = mr/2, integer n. 

Lemma 2 If xiu) is a solution to (H) with x(l) = 7i'/2 and x?/(l) > 1 
then there exists a < < 1 such that xivo) = 0, if Xy{^) < ~1 then there 
exists a < ?/o < 1 such that xivo) = t^- 

A brief outline of the proof: 

We make the substitution y = tan 6 which gives the equation 

cos2^(sin2^Xee + sin2^Xe) = sin2x(^). (28) 

After differentiation of this equation one can get some inequalities on the 
third derivative of Xy these can then be used in order to show that if 
Xe(7r/4) > 2 then xeiO) > Xe(7r/4), < ^ < n/A. ( Note that XeiTr/A = 
2xy(l).) This leads to the existence of a < 6^0 < 7r/4 such that x(^o) = 0. 
A similar argument shows that if xeiT^/A) < —2 then there exists a < < 
7r/4 such that x(^o) = tt. 



-2/ dxxxl = -[cos2x{x)]°^. (26) 
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Lemma 3 If xiv) is a solution to (||) with = tt/2 and < 1 

then < xiy) < v", ?/ > 1. 
An outhne of the proof: 

We denote the known solutions with = 7r/2 by = 7r/2 ± 

(2 arctana; — 7r/2). Using the equation we get if we multiply (^) with Xx and 
integrate from x to 1, we can show that if < 1 then < 

for < a; < 1. Since < x^'i^) < we thus have < x{^) < tt, < a; < 1. 

Proof of the theorem : 

If x{y) is a regular solution to with x(0) = and Xy{^) > then 
the corollary tells us that = tt/2. From lemma 1 together with lemma 
2 it follows that we can not have > 1- Together with the existence 

of the solution xiu) = 2 arctan?/ (which has = 1 and x(oo) = tt) 

we thus conclude that —1 < Xy{^) ^ 1- From lemma 3 it now follows 
that 0<x(?/)<7r,?/>l, this together with Lemma 1 thus tells us that 
< xiy) < ^) for all finite y. A similar reasoning gives — tt < xiv) < if 
Xj^(O) < 0, which completes the proof. 
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Figure Caption 



1 The previously known self-similar solution to the non-linear sigma 
model (solid line) and the first four of the new class of solutions found 
in this paper. Here xiv) is the radial function in the spherically sym- 
metric ansatz, and y — r/tis the self-similarity variable. The values of 
the derivative at the origin are for each solution, Xy — '^i — 21.757, 
Xl = 234.50, xl = 2521.3 and xt = 27102. 
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